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The effects of finite temperature in transport through nanoscopic systems exhibiting uniaxial 
magnetic anisotropy D, such as molecular magnets, adatoms, or quantum dots side-coupled to a 
large spin are analyzed in the Kondo regime. The linear-response conductance is calculated by means 
of the full density-matrix numerical renormalization group method as a function of temperature T, 
magnetic anisotropy D, and exchange coupling J between the molecule's core spin and the orbital 
level. It is shown that such system displays a two-stage Kondo effect as a function of temperature 
and a quantum phase transition as a function of the exchange coupling J. Moreover, additional 
peaks are found in the linear conductance for temperatures of the order of T ~ \J\ and T ~ D. It is 
also shown that the conductance variation with T remarkably depends on the sign of the exchange 
coupling J. 
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I. INTRODUCTION 

Recent progress in experimental techniques that al- 
low for dealing with systems involving single atoms or 
molecules has opened a path for a new generation of 
electronic and spintronic devices Functionality of such 
systems is usually based on their magnetic properties!^! 
In particular, the combination of a large spin and uniax- 
ial magnetic anisotropy makes magnetic adatoms^ and 
single-molecule magnets (SMMs) 8 promising candidates 
as information-storage mediaP^ 

The understanding of transport properties of atomic 
and/or molecular systems exhibiting magnetic anisotropy 
in the whole range of the coupling to electrodes lies 
at the bottom of their potential applications. Espe- 
cially interesting in this context seems to be the limit of 
strong coupling, in which some nontrivial many-body ef- 
fects stemming from the interplay of magnetic anisotropy 
and the Kondo effect are expectedP"2l In particular, it 
turned out that the cooperation of quantum tunneling 
and spin-exchang e proc esses may lead to the pseudo-spin- 
1/2 Kondo effect EEU Moreover, as long as a moderate 
external magnetic field is involved, there are no quali- 
tative differences between the mechanisms of the Kondo 
effect in systems with half- and full-integer spins. 12 It has 
also been suggested that the formation of the Kondo res- 
onance should depend on how the system's total spin is 
modified, i.e. reduced or augmented, upon accepting a 
surplus electron!^ Actually, the Kondo effect can occur 
only in the former case, which corresponds to the antifer- 
romagnetic coupling in the effective spin- 1/2 anisotropic 
Kondo Hamiltonian. Interestingly enough, by changing 
the magnitude of transverse magnetic field one can in- 
duce the oscillations of the Kondo effect, which stem 
from the Berry-phase periodical modulation of the tun- 
nel splitting.^ Finally, when the attached electrodes are 
ferromagnetic, behavior of the tunnel magnetoresistance 



(TMR) in the Kondo regime for systems under discussion 
is expected to be significantly different 13 * 14 ! from that for 
typical magnetically isotropic quantum dotsP^E^ 

Although the observation of the Kondo-related fea- 
tures in systems displaying magnetic anisotropy is ex- 
perimentally challenging as it requires cooling the system 
down to very low temperatures, several successful mea- 
surements have been recently reported {EMU Generally, 
variation of temperature in a nanoscopic system reveal- 
ing the Kondo correlations results in a dramatic change 
of its transport properties. In quantum dots, lowering 
the temperature T below a characteristic energy scale - 
the Kondo temperature Tk - is accompanied by an in- 
crease of the conductance to its maximum value, which, 
on the other hand, becomes gradually suppressed with in- 
creasing TP^^In the Kondo regime, the dependence of 
the linear conductance on temperature is then a univer- 
sal function of T/Tk^^ Moreover, it was shown very 
recently that the conductance of quantum dots coupled 
to ferromagnetic leads or subject to an external magnetic 
field also exhibits univer sal fe atures with respect to the 
(effective) magnetic field pSE3 

The temperature dependence of the Kondo effect be- 
comes more complex when the system consists of more 
impurity spins, and the competition between the spin- 
exchange processes due to tunneling of electrons and 
the interaction between the constituent spins is possi- 
ble. Even in the conceptually simplest case involving two 
exchange coupled spin-1/2 impurities, e.g. as in quan- 
tum dots containing an even number of electrons,* 28 * 29 ! 
there are several different scenarios regarding the inter- 
impurity exchange interaction j ! 3Q " 35 l In the low tem- 
perature limit, the impurities form a singlet (5 = 0) 
for large antiferromagnetic J and the effect of conduc- 
tion electrons on the system is weak, while the high-spin 
triplet (5 = 1) ground state develops for the ferromag- 
netic coupling, which then can be screened. A similar 
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scenario is also relevant for side-coupled double quantum 
dots in single spin regime, with only one dot coupled 
directly to conduction electrons. In such a setup, de- 
pending on the strength and sign of the spin exchange 
interaction between the two dots, a two-stage Kondo ef- 
fect and Kosterlitz-Thouless quantum phase transition 
can occurP^ESl 

As the magnitude of the impurity spin becomes larger 
than 1/2, a new energy scale related to magnetic 
anisotropy enters the problem. In the case of a single 
anisotropic Kondo impurity with spin S > 1/2 coupled 
to a single conduction-electron channel, which accurately 
represents the situation o f a m agnetic adatom deposited 
on a nonmagnetic surface,^ * 19 * 20 ! it has been shown that 
at low temperatures the s pin is ultimately always subject 
to complete compensationP^^In this paper we consider 
more complex systems exhibiting an uniaxial magnetic 
anisotropy, e.g. adatoms, quantum dots side-coupled to 
a large spin 5, as well as molecular magnets, where ad- 
ditionally the charge state of the system can be changed 
owing to electron tunneling processes. The central aim 
is then to study the finite-temperature transport proper- 
ties of such systems, focusing on the Kondo regime. Un- 
derstanding the behavior of system at different temper- 
atures is of great importance, as many experiments are 
actually carried out in the cross-over regime, T ~ Tk, 
when neither Fermi liquid (T <C Tk) nor perturbative 
(T ^> Tk) descriptions are applicable. Our analysis is 
based on the full density-matrix numerical renormaliza- 
tion group (NRG) method ] 43 " 45 1 which allows for calculat- 
ing transport at any temperature in an essentially exact 
way. 

II. THEORETICAL DESCRIPTION 

The key features of systems under consideration - such 
as magnetic adatoms, quantum dots coupled to localized 
magnetic impurities, and SMMs - can be reproduced by 
a model consisting of a single conducting orbital level 
(OL) exchange-coupled with strength J to a magnetic 
core (spin S) subject to uniaxial magnetic anisotropy 
D, see Fig. [T] This system will be referred to as mag- 
netic quantum dot (MQD) and its Hamiltonian has the 

fori jmn3] 

Hmqd = Hol - DS 2 Z -Js-S. (1) 

The orbital level is described by Hol = £ J2 a n a +U n^n^, 
where n a = cj.c a is the occupation operator, with cl(c a ) 
creating (annihilating) a spin-a electron of energy e in 
the orbital level, while U accounts for the Coulomb en- 
ergy of two electrons of opposite spins residing in the 
orbital. Furthermore, the second term of the Hamil- 
tonian ([TJ characterizes the magnetic anisotropy of the 
core, where D denotes the uniaxial magnetic anisotropy 
constant, and S z is the zth component of the MQD's in- 
ternal spin operator S. The present discussion is limited 
only to the case of systems exhibiting magnetic bistability 




Magnetic Quantum Dot 



FIG. 1. (Color online) Schematic depiction of the system 
under consideration. It consists of a single orbital level (OL) 
tunnel-coupled to external leads, with coupling strengths 
Yl and Yr for the left and right leads, and additionally 
exchanged-coupled to a spin S, with J denoting the strength 
of exchange interaction. 

(D > 0). Finally, the exchange interaction between the 
magnetic core of a MQD and the spin of an electron oc- 
cupying the orbital level, given by s = \ ^ a(J , c\a aa ,c a , , 
with a = (cr x , a y , cr z ) standing for the Pauli spin opera- 
tor, is expressed by the last term of Hmqd- In general, 
the J-coupling can be either of ferromagnetic (FM for 
J > 0) or antiferromagnetic (AFM for J < 0) type. 

The magnetic quantum dot is assumed to be tunnel- 
coupled to two identical electrodes only via the or- 
bital level, see Fig. [I] and electrons in the qth elec- 
trode [q = (L)eft, (i?)ight] are modelled by, = 

Xlcr I^w ^ e ea Jcr( e ) a gcr(e), with a qa (e) being the relevant 
creation operator, and W the band half-width. The 
electron tunneling processes between the MQD and elec- 
trodes are described by 

fttun = y\/— / de[al a (e)c a + cla qa (e)], (2) 

qa V ^ J W 

where Y q represents the strength of coupling of the orbital 
level to the qth lead. 

Conceptually, the model considered is equivalent to a 
single-level quantum dot which - if occupied by a single 
electron - is exchanged coupled to a large-spin magnetic 
impurity subject to magnetic anisotropy! 3 -^ Moreover, to 
some extent it can be regarded as an alternative to a two- 
impurity Kondo model where only one impurity couples 
directly to conduction band, or to a double quantum dot 
system in a T-shape geometry! 37 " 40 ! From this point of 
view our model will also exhibit a two-stage Kondo effect 
and a quantum phase transition, as shall be discussed in 
next section. 

Generally, the physics of the Kondo effect is essen- 
tially determined by the number of conduction-electron 
channels to which the system is coupled, as in order to 
completely screen a spin S one needs 2S channels. 4 ^ In 
this context, high- spin molecular devices are unique as 
they commonly operate in the regime where effectively 
only one channel plays a role! 21 * 35 ] More precisely, even 
if the device is coupled to multiple leads, each of such 
junction usually supports only a single conduction chan- 
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nel, and additionally the couplings are typically charac- 
terized by a strong asymmetry. As a result, the rele- 
vant Kondo energy scale is determined by the strongest 
coupling, since the Kondo temperatures associated with 
weaker couplings are exponentially small, =^ and hence 
negligible under typical experimental conditions. In the 
model considered here the situation is even simpler, since 
the orbital level of magnetic quantum dot couples only 
to an even linear combination of the electron operators 
in the left and right leads, with a new coupling strength 
T = T^+r^, while the odd combination is completely de- 
coupled. This basically means that we can limit our dis- 
cussion to the case of a single conduction-electron chan- 
nel. 

In the following, we will study the linear-response 
transport properties of magnetic quantum dots on var- 
ious parameters of the system. The main quantity we 
are interested in is the temperature-dependent linear 
conductance G, which we ca lculate from the Landauer- 
Wingreen-Meir formula f 48 * 49 * 
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where f(uo) denotes the Fermi-Dirac distribution func- 
tion, while A a (uj) is the spin-dependent spectral func- 
tion of the orbital level, and 2T L T R /(T L + T R ) = T/2 
in the present situation. The problem of determining the 
MQD's transport features corresponds then essentially to 
finding the spin-resolved spectral function A a (ui), which 
in the present work is obtained by means of the Wilson 's 
numerical renormalization group (NRG) method.E^Mlln 
particular, we employ the recent idea of a full density ma- 
trix,^ which allows for reliable calculation of static and 
dynamic properties of the system at arbitrary temper- 
atures.^! For the present problem, to obtain decent re- 
sults, the Z7 c harge(l) x ^spin(l) symmetry was exploited, 
the discretization parameter A = 1.8 was used and we 
kept Nk = 1200 states during calculations. 



III. RESULTS AND DISCUSSION 

In this paper we consider the transport features of a 
prototypical magnetic quantum dot characterized by the 
spin S = 2. The other parameters are typical of molec- 
ular systems, see the caption of Fig. [2] In order to dis- 
cuss the influence of finite temperature on the Kondo ef- 
fect, we introduce the Kondo temperature Tk, to which 
other parameters will be compared whenever it is useful. 
The Kondo temperature is defined here as the half-width 
at the half-maximum of the normalized linear conduc- 
tance G/Gt,j,d=o as a function of T, where Gt,j,d=o 
is the conductance for T = D = J = 0, which yields 
T K /U ^5-10- 3 . 
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FIG. 2. (Color online) The dependence of the linear con- 
ductance G on the energy of the orbital level (OL) e and 
the temperature T in the case of (a,c) ferromagnetic (FM) 
and (b,d) antiferromagnetic (AFM) exchange coupling J for 
\J\/Tk = 3 and D = 0. The variable Q represents the average 
number of electrons that occupy the OL. Here, Gt,j,d=o cor- 
responds to the conductance calculated for T = D = J = 0. 
Parts (c)-(d) present the relevant cross-sections of the density 
plots at the electron- hole symmetry point (e = —U/2) for dif- 
ferent values of the J-coupling, as indicated. The parameters 
of the system are as follows: U = 10 meV, T/U = 0.1 and 
T K /U ^5-10~ 3 . 



A. The case of zero magnetic anisotropy (D = 0) 

Before the influence of the magnetic anisotropy on the 
temperature dependence of the Kondo effect is analyzed, 
it is instructive to consider the situation with vanish- 
ing magnetic anisotropy, D = 0. The corresponding re- 
sults are shown in Fig. [2j which presents the dependence 
of the linear conductance on orbital level position and 
temperature for different exchange couplings J. It can 



4 



be seen that when the orbital level is occupied by even 
number of electrons, which corresponds to an empty or 
fully occupied level, the conductance is generally sup- 
pressed and determined only by elastic cotunneling pro- 
cesses, see Figs. [2|a,b). On the other hand, for an odd 
occupation of the orbital level (single spin regime), the 
Kondo effect should in general develop at low tempera- 
tures. Now, however, the exchange interaction with the 
core spin S comes into play. In principle, two distinctive 
cases with respect to the J-coupling sign can be recog- 
nized. For the ferromagnetic exchange coupling, the low 
temperature behavior of the system is governed by the 
competition between the coupling J and hybridization 
with electrodes T. The former interaction tends to stabi- 
lize the high- spin state, whereas the later one leads to the 
screening of the orbital level's spin. It has been shown 
for a system including two spin- 1/2 impurities, that the 
Kondo effect dominates even if the exchange coupling is 
significantly larger than the hybridization ! 36 ! 51 ! This can 
be also seen in Fig. |5|c), where irrespective of the mag- 
nitude of J, G tends to its maximum value in the zero 
temperature limit. Nevertheless, the exchange interac- 
tion still plays a prominent role, because it is responsible 
for the reduction of the Kondo temperature, as compared 
to a simple spin- 1/2 magnetic impurity system. 



The situation is qualitatively different for antiferro- 
magnetic J-coupling. Investigations of high-spin two- 
impurity (with spins Si and S2) Kondo models charac- 
terized by the antiferromagnetic inter-impurity exchange 
interaction and only one impurity (Si) directly coupled 
to a conduction band have revealed a two-stage Kondo ef- 
fect which is a generic feature for all S2 ^ Si models.^ In 
particular, a two-stage Kondo process should take place 
for the inter-impurity exchange coupling smaller than the 
energy scale Tk associated with the Kondo screening of 
the directly coupled impurity. For Si = 1/2, the temper- 
ature T K at which S2 becomes screened is then exponen- 
tially smaller, as the process occurs due to interaction 
of Si with a Fermi sea arising as a consequence of the 
first screening staged The two-stage Kondo effect can 
be nicely seen in Fig. [2^d), which shows the tempera- 
ture dependence of linear conductance for different an- 
tiferromagnetic exchange coupling. Let us have a closer 
look at the curve corresponding to \J\/Tk — 0.2. For 
T < Tk, the conductance starts increasing due to the 
Kondo effect associated with the screening of the orbital 
level spin by conduction electrons. However, at suffi- 
ciently low temperatures the energy scale related with the 
second stage of screening becomes relevant and the con- 
ductance becomes suppressed. Moreover, when varying 
the strength of the exchange interaction at zero temper- 
ature, the system undergoes a quantum phase transition. 
This phase transition is similar to a singlet-triplet tran- 
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sition observed in multilevel quantum dotsj 32 ' 33 1 35 1 52 ^ 
and for finite temperature, magnetic field or anisotropy 
turns into a crossover. 



FIG. 3. (Color online) The normalized linear conductance 
G as a function of temperature T and exchange coupling pa- 
rameter J in the case of (a,c) ferromagnetic (FM) and (b,d) 
antiferromagnetic (AFM) type of the coupling between mag- 
netic core and orbital level shown in the presence of uniaxial 
magnetic anisotropy, D/U = 10 -4 (D/Tk = 0.2). Curves 
in (c)-(d) represent cross-sections of figures (a) and (b), re- 
spectively, for indicated values of |J|, but now each curve is 
normalized to the corresponding Gt=o instead of Gt,j,d=o as 
in (a)-(b). The other parameters are the same as in Fig. [2] 
with e -U/2. 



B. The case of finite magnetic anisotropy (D > 0) 

In the case of finite magnetic anisotropy D, the situ- 
ation becomes more complex since now the degeneracy 
of spin multiplets is partially lifted. More precisely, in 
the absence of magnetic anisotropy D, the ground state 
of the MQD is the spin multiplet S + \ (S — \) for fer- 
romagnetic (antiferromagnetic) exchange interaction J. 
In the case of finite magnetic anisotropy (D > 0), the 
ground state becomes two-fold degenerate and consists 
of the lowest-weight and highest-weight components of 
the above multiplets, respectively.^ 

Figure [3] presents the dependence of the system's lin- 
ear conductance G on the exchange coupling J and tem- 
perature T, for the ferromagnetic and antiferromagnetic 
types of the J-coupling. As long as \J\ < Tk, the tem- 
perature behavior of the conductance does not diff er from 
that for a typical single- level quantum dot) 18 | 27 | ^ see the 
parts (a) and (b) in Fig. [3] The conductance also does 
not depend on the sign of the parameter J, compare the 
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FIG. 4. (Color online) The normalized linear conductance G as a function of temperature T and uniaxial magnetic anisotropy 
D for \J\/Tk = 0.4 in (a)-(d) and \ J\/Tk = 2 in (e)-(h). The other parameters are the same as in Fig. [5] with e — —U/2. 
Dashed lines in (a)-(b) and (e)-(f) represent maximal value of the uniaxial magnetic anisotropy constant D typical to molecular 
magnets, for further details see the text. 



solid lines (for \ J\/T K = 0.1) in Fig. gc)-(d). In such 
a limit, electrons tunneling through the orbital level are 
hardly affected by the presence of the MQD's magnetic 
core. However, as \J\ > Tk, the screening of the orbital 
level's spin by conduction electrons is suppressed due to 
the strong exchange interaction of the orbital level's spin 
with the magnetic core. One can then observe that while 
for T « Tk the conductance G just decreases mono- 
tonically with increasing | J|, some additional features of 
G emerge for T ~ Tk, see Figs. [3]^c)-(d). This stems 
from the fact that for large \J\ the influence of the mag- 
netic core on the orbital level cannot be neglected, as the 
tunneling processes actually take place via the molecular 
spin states formed due to the J-coupling. 

Since transport properties of the magnetic quantum 
dot in the linear response regime depend basically on the 
system's ground state, thus for \J\ > Tk, it is the sign 
of the exchange parameter J that determines the sys- 
tem's spin multiplets that plays the dominant role, i.e. 
5+1/2 for ferromagnetic J and S — 1/2 for antiferro- 
magnetic J. Specifically, at low temperatures the MQD 
with ferromagnetic J is found to occupy the doublet state 
S z = ±5/2 of the spin multiplet 5 = 5/2, whereas under 
the same conditions the system with antiferromagnetic J 
prefers the state S z = ±3/2 belonging to the spin mul- 
tiplet S = 3/2. Noting this, the characteristic behavior 



of the linear conductance for temperatures around Tk in 
Fig. [3] can be explained straightforwardly. In general, the 
conductance displays additional features (peaks) when- 
ever the number of MQD's states participating in elec- 
tronic transport changes due to increased temperature. 
For this reason, as soon as T ~ D, the neighboring states 
of the same spin multiplet become accessible for electron 
tunneling processes, and when T reaches the value of | J|, 
T ~ |J|, also the states of the other spin multiplet en- 
ter into consideration. At these temperatures the linear 
conductance displays additional peaks, see Figs. [3] (c)- 
(d). Nevertheless, it should be emphasized that energies 
of the MQD's states with the orbital level occupied by a 
single electron generally depend both on J and 

Additional feature visible in Fig.|3]is an asymmetry be- 
tween the ferromagnetic and antiferromagnetic coupling 
J. Note, however, that the value of Gt=o> to which the 
conductance in Figs. [3] (c)-(d) is normalized, is different 
in both cases, and the difference is visible only if | J\ > T, 
otherwise the thermal fluctuations smear the difference 
between the cases of ferromagnetic and antiferromagnetic 
exchange coupling, Figs. [3] (a)- (b). This is related to dif- 
ferent spin multiplets relevant for spin flip processes that 
drive the Kondo effect at low temperatures.^ 

Let us now focus on the effects due to the magnetic 
anisotropy T>, see Fig. |1J It is worth noting that the 
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FIG. 5. (Color online) The normalized linear conductance 
G as a function of temperature T and orbital level position e. 
Parts (c)-(d) present the relevant cross-sections of the density 
plots for different temperatures, as indicated. The parameters 
are the same as in Fig. [2] with \ J\/T K 3 and D/T K 0.2. 



anisotropy constant in SMMs can take different values, 
e.g. in the case of a Fe4 molecule one finds D/U ~ 
5 • 10 -3 . This is actually one of the highest values of 
D observed in SMMs. On the other hand, in magnetic 
adatoms, such as e.g. Fej^this constant can be as large 
as D/U ~ 0.1. Therefore, we study the effects resulting 
from magnetic anisotropy and finite temperature for a 
broad range of anisotropy constant. Figure [4] shows the 
linear conductance vs. temperature and anisotropy con- 
stant for both positive and negative exchange coupling 
J. It can be seen that significant differences between 
the case of the ferromagnetic and antiferromagnetic J- 
coupling occur for D <C |J|. When the anisotropy is 
weak, the behavior of linear conductance is the same 
as discussed previously in the case of D = 0, i.e. at 
low temperatures the system is in the underscreened 
Kondo regime for ferromagnetic J-coupling, while for the 
antiferromagnetic J-coupling one observes a two-stage 
Kondo effect. When the anisotropy increases, see e.g. 
the case of D = Tk, both the underscreened and two- 
stage Kondo effects become generally suppressed. More- 
over, the temperature dependence of the normalized lin- 
ear conductance displays a qualitatively different behav- 
ior for the ferromagnetic and antiferromagnetic type of 
the J-coupling. While in the case of ferromagnetic J 
the increase of T generates a drop of the conductance, 



see Figs. [4ja,c), the opposite trend appears for antifer- 
romagnetic J, see Figs. Qb,d). This is mainly related 
with the fact that for a given value of |J|, the zero- 
temperature conductance is much smaller for the antifer- 
romagnetic case than for the ferromagnetic one. Inter- 
estingly enough, in both situations one can notice a peak 
at T ~ I J|. For this temperature, the molecular states of 
another multiplet corresponding to the singly occupied 
orbital level become available for transport. Because of 
it, the overall rate of spin- flip processes is increased due 
to thermal fluctuations and an additional resonance in 
the conductance appears. Finally, when D > |J|, the 
differences between J < and J > actually disappear, 
see the dotted curves for D/T K = 10 in Figs, gc)-(d), 
which are almost identical. Note that the influence of 
magnetic anisotropy of the core spin on the transport 
properties of the system is only present for considerable 
exchange couplings | J| > Tk- 

It is also interesting to analyze the temperature depen- 
dence of the linear conductance in the presence of mag- 
netic anisotropy and for varying position of the orbital 
level £, as shown in Fig. [5] The orbital level position can 
be experimentally tuned by a gate voltage. By changing 
£, the orbital level becomes consecutively occupied with 
electrons, see Fig. [5] where Q denotes the average number 
of electrons. At low temperatures, T <C T, transport for 
even occupations is mediated by cotunneling processes 
and conductance is low. For the odd occupation (Q = 1) 
and T < Tk, the Kondo effect should generally develop. 
It is however suppressed due to finite exchange coupling 
and magnetic anisotropy. It is very instructive to com- 
pare Figs. [5ja) -(b) with Figs.[2ja)-(b), where actually the 
same dependence is plotted for finite anisotropy (Fig. [5| 
and for D = (Fig. 2|. It can be seen that the main 
difference occurs in the case of ferromagnetic exchange 
coupling, when finite anisotropy leads to the suppression 
of the Kondo effect and conductance is much lower than 
in the case of D = 0. In addition, when T ~ |J|, a 
resonance due to thermally-activated spin-flip processes 
through other spin multiplets occurs. This effect is more 
visible in the case of antiferromagnetic J. On the other 
hand, once T > |J|, the conductance starts decreasing 
and the difference between the cases of positive and neg- 
ative J is diminished. 

Finally, we study the normalized linear conductance 
as a function of the anisotropy constant D and exchange 
coupling J for a constant temperature, T /Tk = 10 -2 , see 
Figs. [6ja)-(b). One can see that the effect of exchange 
interaction becomes visible when |J| > Tk- If this is 
the case, then the transport properties also depend on 
the anisotropy constant D. For a fixed value of | J|, with 
|J| ^ Tk, the conductance becomes suppressed by in- 
creasing the anisotropy constant in the case of ferromag- 
netic exchange interaction, while for antiferromagnetic 
J, the behavior is just opposite - there is an increase of 
G. This tendency is also nicely seen in the cross-sections 
shown in Figs. [6]^c)-(d) for a few anisotropy constants 
D. Generally, the conductance drops with increasing J 
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FIG. 6. (Color online) Dependence of the normalized lin- 
ear conductance G on the magnetic anisotropy D and the 
exchange coupling J of (a) the ferromagnetic (FM) and (b) 
antiferromagnetic (AFM) type for T/T K = 0.01. Parts (c) 
and (d) show the dependence of G on (c) ferromagnetic and 

(d) antiferromagnetic exchange coupling J for T/Tk — 10 -2 
and for a few anisotropy constants D, as indicated. Parts 

(e) and (f) show the same calculated for T —> 0. The other 
parameters are the same as in Fig. [2] with e = -U/2. 



irrespective of the type of exchange interaction. This is, 
however, because the curves present G vs. J at finite 
T/Tk = 10 -2 . In the case of antiferromagnetic exchange 
coupling there is a quantum phase transition as J is var- 
ied. This can be seen in Fig. |6jf), which was calculated 
for temperature T — >• 0. This transition turns into a 
cross-over in the case of finite T and D, see Figs. [6^d) 
and (f). In the case of ferromagnetic exchange interac- 
tion in the absence of magnetic anisotropy and for T — >• 0, 
the conductance does not depend on J and equals 2e 2 //i, 
see Fig. J^e). For finite anisotropy D, the degeneracy of 
the ground state multiplet S + \ is lifted and the Kondo 



effect is suppressed once \J\ > Tk- On the other hand, 
if the temperature is finite, then at certain value of fer- 
romagnetic exchange interaction the conductance drops 
since the screening temperature depends on J and the 
condition T < Tk is not met any more. 



IV. SUMMARY 

We have studied transport properties of magnetic 
quantum dots coupled to external leads in the Kondo 
regime. The analyzed system consisted of a spin S = 2 
exchange coupled to a single orbital level that was di- 
rectly tunnel-coupled to electrodes. In particular, we 
have focused on the dependence of the linear-response 
conductance of the system on temperature, orbital level 
position, magnetic anisotropy and exchange coupling. 
The calculations were performed with the aid of full 
density-matrix numerical renormalization group method. 
In the absence of magnetic anisotropy, the model studied 
generally exhibits a two-stage Kondo effect and an un- 
derscreened Kondo phenomenon, depending on the sign 
of the exchange coupling. We have shown that these 
two effects become generally suppressed in the presence 
of magnetic anisotropy, if the exchange coupling is suf- 
ficiently strong, \J\ > Tk- We have also shown that 
the temperature dependence of the conductance depends 
on the type of exchange coupling J. For ferromagnetic 
coupling, the linear conductance was found to generally 
decrease with T, while for antiferromagnetic coupling, 
the conductance displayed a maximum for T ~ |J|. In 
addition, the conductance variation with the exchange 
coupling reveals a quantum phase transition for antifer- 
romagnetic J-coupling, which turns into a crossover in 
the case of finite temperature and magnetic anisotropy. 

The model studied in the present paper can be used 
to describe transport through single-molecule magnets, 
adatoms or quantum dots exchange-coupled to a large 
spin S. Our results and predictions may be thus relevant 
for a wide class of molecular devices. Finally, we note 
that molecular devices are more suitable for observing 
Kondo-related effects, since they provide larger energy 
scales, which translates into higher Kondo temperatures, 
more easily achievable in experiments.^ 
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